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Dissipation statistics of a passive scalar in a 
multidimensional smooth flow * 

Andrea Gamba^ and Igor V. Kolokolov^ 


Abstract 

We compute analytically the probability distribution function V (e) 
of the dissipation field e = (V0)^ of a passive scalar 9 advected by a 
d-dimensional random flow, in the limit of large Peclet and Prandtl 
numbers (Batchelor-Kraichnan regime). The tail of the distribution 
is a stretched exponential: for e —> oo, In'P(e) 
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1 Introduction 

Intermittency and strong non-gaussianity of developed turbulence are most 
clearly reflected in the peculiar structure of the observed probability distri¬ 
bution functions (p.d.f.) of the gradients of the turbulent held. A typical 
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logarithmic plot of the gradients p.d.f. is concave rather than convex, show¬ 
ing a strong central peak and slowly decaying tails. Rare strong fluctuations 
are responsible for the tails, while large quiet regions are related to the cen¬ 
tral peak. In particular, such p.d.f. were observed for the square gradients 
£ = (V6*)^ (dissipation held) of a scalar held 9 passively advected by an 
incompressible turbulent how. 

From the theoretical standpoint, most rigorous results in the theory of 
turbulent mixing have been obtained so far within the framework of the 
Kraichnan model B i, describing the advection of a passive scalar by a 
random velocity held, delta-correlated in time. The exact solvability of the 
Kraichnan model suggests that it can play in turbulence a role similar to 
that played by the Ising model in the theory of critical phenomena [|. 

In particular, the p.d.f. of the dissipation held £ was recently calculated 
analytically in the one-dimensional and two-dimensional ^ cases. In 
this paper we extend the result of Ref. to arbitrary space dimensions 
d. Our techniques are based on a combination of functional integration and 
group-theoretical methods, hrst introduced by one of us 0 in the context of 
quantum magnetism. These techniques allowed hrst to compute exactly the 
p.d.f. of the passive scalar in the two-dimensional, non-dissipative case p. 
The multidimensional case was considered in P, |T^. In the dissipative case. 


an essential ingredient is the time-separation method introduced in in a 
one-dimensional context, and afterwards exploited in Ref. |^. The point is 
that in developed turbulence there appears a natural small parameter, 1 / Pe, 
where Pe is the Peclet number, measuring the relative strength of advection 
with respect to dihusion. One would like to use this parameter to develop an 
asymptotic theory for quantities like the p.d.f. V{e) of the dissipation held 
£. However, such quantities are essentially non perturbative in 1/Pe. For 
instance, £ is exactly zero without dihusion, but it has a non-zero limit as 
1/Pe 7^ 0, Pe —*• cxo. ^ 


^The similarity should be understood in a loose sense. The Kraichnan model can 
be viewed in itself as a sort of mean-field approximation to a full theory of developed 
turbulence @]. 

^It is worth noting that without dissipation stationary limits for the p.d.f.’s of the 
passive scalar field, of its gradients, etc., do not exist. In the case of the direct cascade 
considered in this paper an arbitrary small but finite dissipation constant k provides 
stationary distributions at large times for all the physical quantities. However, there are 
cases corresponding to the dynamical inverse cascade when dissipation does not lead to 
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Following Refs. Bi. we show that an appropriate dynamical formalism 
natnrally introdnces two time-scales, a short one related to stretching, and 
a long one related to diffusion. The time scale of the stretching fluctuations 
is of the order of the inverse of the maximum Lyapunov exponent, while the 
whole time of stretching is InPe times larger. 

The paper is organised as follows: 

In Sec. 1^ we dehne the problem, recall the statistical and kinematical 
concepts involved, introduce the functional transformation which makes the 
problem solvable and write down the corresponding functional Jacobian. 

In Sec. ^ we expose the time-separation method, showing that the com¬ 
putation of 'P{e) can be split in two parts, which we will denote as small- 
and large-times averaging. Suggestively, this procedure can be thought of as 
a two-step application of a renormalisation group transformation. 

In Secs. ^ and we perform respectively large- and small-times averag¬ 
ing, reducing the computation of the corresponding functional integrals to 
solvable auxiliary quantum-mechanical problems. 

In Sec. 1^ we put together all the pieces of the computation and present the 
resulting p.d.f. V{e), whose asymptotic behaviour is discussed as a function 
of the space dimension d. 

2 Statement of the problem 

The advection of a passive scalar held 6{t, r) by the incompressible smooth 
how v(t, r) in d-dimensional space is governed by the transport equation 

{dt + v-W - KA)e = (p, (1) 

where r) is an external source, k is the dihusivity and we assume a null 
initial condition at t = —oo. We assume that both v(t, r) and (p{t, r) are 
Gaussian independent random functions, J-correlated in time [|^ |^, |^, Q 
so that, hrstly: 


(0(ti, ri)0(t2, ra)) = J(ti - t2)x(ri2), (2) 

where ria = |ri — ra], the correlation ^(ria) decays on the scale L, and the 
constant Pa = x(0) is the production rate of 6"^. We consider large Prandtl 

stationary p.d.f.’s for some observables [pi]. 
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numbers, which correspond to a large viscosity-to-diffusivity ratio, so that in 
the viscous interval 


,(ti, ri)vi3(t2, r2)) = S(ti- t2) VoSap - - Targ) 


(d- 1)D 
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( 3 ) 

Isotropy of the velocity statistics is here assumed. The representation (^) is 
valid for scales smaller than the velocity infrared cut-off (which is supposed 
to be the largest scale of the problem), since represents the two hrst terms 
of the expansion of the velocity correlation function in powers of r/Lu, and 
2) ~ Vo/L^. We assume also that the inequality Pe^ = 'DL‘^/2k 3> 1 holds, 
guaranteeing the existence of a convective interval of scales rdig r L 
where the dominant effect to be taken into account is the stretching of the 
velocity held. Here rjiff = 2y/ k/{ d — 1)2) is the mean diffusion length. 

It follows from (|]) that the correlation functions of the traceless strain 
held, dehned as aap = / 3 Va, are r-independent: 


(*^0/3(H(^ 2 )) 2D [(d -|- l)dQ:^d/3i/ d(tl ^ 2 ) • (4) 


This means that the strain held aa/s can be treated as a random function of 
time t only. To exploit this property it is convenient to pass to a comov¬ 
ing reference frame, that is, to a coordinate frame whose origin follows the 
motion of a Lagrangian particle of the huid, so that the velocity held can 
be approximated by v ~ d'(t)r, with {aa/sit)) = 0 0. With this substitu¬ 
tion, (|ID can be solved by Fourier transforming in the space coordinates and 
integrating along characteristics: 


dk(t) 



dt' exp 


nk- f dtiW{t,ti)W^{t,ti)k 

Jt' 


( 5 ) 


where 

W(t,F) = Texp 
is the solution of the evolutive problem 



( 6 ) 


Fm = ff(()VV((). W((') = l, (7) 


and T is the chronological ordering operator. We will study stationary sta¬ 
tistical properties of the passive scalar at the hxed time moment t. It is 
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convenient to use the symmetries of the measure of averaging over a(r) with 
respect to the transformations (T(r) —<5 ‘(—t) and (T(r) ^ —d(r). Chang¬ 
ing now variables in the time integral in (^) we arrive to the expression: 


^k(^) = / dt',W ^’"’"(C)k) exp —Kh drVV ^(r)>V . (8) 


Here VV(t) = W{t, 0) and invariance of the random force 0 averaging with 
respect to time inversion was used. Let us now consider the statistics of the 
dissipation held 

£ = 

The p.d.f. of the e held can be written as 

1 p+ioo p Jd™ 

= (9) 

with 

Jo (k->V(t)m)' exp (^-_k-A(t)kJ (10) 

and ^ 

A(t) = W^t) >V-^(r)W-i’T(r)dr^ W^(t) (11) 

Here the integration over the auxiliary vector parameter m takes care 
of combinatorics and summation over vector indices. The moments of the 
dissipation held £ can be recovered as 

/ d'^777 

( 12 ) 


The average (.. .)„ over the traceless random strain held a has to be per¬ 
formed with the probability measure T>d{t) exp{—Sa) corresponding to (^: 


^H-OO 




2d{d + 2)T) Jo 


[{d + l)Tr((j(j'^) -|- Tr((j^)]dt. 


(13) 


The matrix A{t) is invariant under right local (time dependent) rotations 
VV(t) —> VV(t)7^(t), and transforms covariantly under left local rotations 
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W(t) ^ TZ{t)W{t). The quantity Q is invariant under left local rotations, 
which for any fixed t can be absorbed in the d'^k integration. The action 
S'ct is invariant under both left and right global rotations and under the 
transformation d ^ d'^. Together with isotropy of the pumping function 
X, invariance of and A under right rotations implies that Q = Q(m) is a 
function of the square modulus 2 ; = m^, so that we can substitute m —»• ^/zno 
in (13), with 

/ 1 \ 

0 


no = 


(14) 


VO/ 


The m-integration measure in (|[) and (^) can then be substituted by 

1 


2:2 "'d 2 ; 


(15) 


nd/2) 

After these substitutions Q maintains invariance under left local rotations of 

W. 

Let us perform the Iwasawa decomposition of the evolution matrix TV: 

W{t) = RDT-^, (16) 

where, for any fixed t, R is an SO{d) rotation matrix, D a diagonal matrix 
with detD = 1, and T an upper triangular matrix with I’s on the diagonal. 
The action takes then the form: 

1 r°° f k ^-I 

= 


1 roo , . _i 1 

df|Tr(DD )^ +-Tr[D2(t-iT)D-2(t-iT)^] 


2dD 


d 


d + 2 


Tr[R-^R- (DT-^TD 


—iTA—I n j2 


(17) 


where (.. .)j^ denotes the antisymmetric part. 

The invariance of Q under local rotations of W implies that Q does not 
depend on R, so that the R variables can be integrated out from the very 
beginning. The resulting effective action is equal to Sa with the last trace- 
term dropped. It is then convenient to parametrize 


D 


/ eP\ 0, 


jPi _ 

0 / 


opl 


0 \ 
0 


V 0, 0, 


• • • 5 


. . . , 


g“(PlH- ^Pd-l) 
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and 


/ 

1, 

hi, 

h2, • 

hd-1 \ 



0, 

1, 

^23, . 

^2,d 

, (19) 

V 

0,’ 

0,’ 

0,’ . 

1 / 



denoting by r]j = the elements of the hrst row of T, which, together 

with the pj, will be seen to be the only relevant dynamical variables. In order 
to have more compact expressions it will be also convenient in what follows 
to denote the snm pi + ■ ■ ■ + pd-i by —pd and to set Xjj = 1. 

The initial condition >V(0) = 1 implies 


p,(0) = 0, P,(0) = 0, J = (20) 

Xmn{^) = 0, 2 < m < n < d. (21) 

With the snbstitntion a (R, D, T) we introdnced in the fnnctional integral 
a set of new “collective” coordinates, by means of a non-linear, non-local 
variable transformation. The Jacobian of the transformation is given by 


J = Jn\- J\ = W exp 
t=o 


■ d-l 

j)pj{t) 

- 


• exp 



d-l 

'^{d- j)pj{T)dT 

i=i 


( 22 ) 

The details of the derivation are given in the Appendix. The local part Ji of 
the Jacobian coincides with the one compnted in Ref. |^. 

With the decomposition (0)> and after the integration of the R variables, 
the probability measnre takes the form 


d-l 


n VXmnJip) exp{-Ses[p,r],X]) , 


2<m<n<d 


with 


S-eff = 


/*CXD 

2^ 


j=l ^ \<i<i<d 


and ^ = T ^T. Let ns also rewrite (pR[) as 


D 


whereby hxing the residnal left-rotation symmetry. 


(23) 


k-A(t)k , (24) 
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3 Time separation 


We shall now expose a method that allows to compute V{e) in the limit 
of large Peclet numbers, he., in a regime characterized by a large ratio of 
advection to diffusion at the pumping scale. 

The moments of the dissipation held e have hnite, non zero limits as 
diffusivity tends to zero. Thus, V{e) is a non-perturbative quantity in the 
small parameter 1/Pe. 

The time separation method, which was introduced in Ref. in a one¬ 
dimensional context and generalized to d = 2 in [^, is the proper tool for 
performing such a non-perturbative calculation. As a matter of fact, the 
dynamical formalism itself reveals the presence of two different time scales: a 
short one, related to diffusion, and a long one, related to advective stretching. 
Taking properly into account the two time scales it is possible to exploit non- 
trivially the presence of the large parameter Pe and to develop an asymptotic 
theory that captures the dominant term in P{e) with respect to Pe. 

The starting point is that, as is clear from the structure of (|^) , the 
main contribution to the moments is obtained in the limit Pe ^ oo for 
Pi ~ InPe. On the other hand, one has the following result (see 0) for the 
hierarchy of the Lyapunov exponents: 

(Pi(^))<J = ^^id- 2j + l)t = Xjt , (25) 


showing that the relevant values of pi are realized for times 

2 


t ~ 


In Pe. 


d(d-l)D.. 

For large times the elements of the matrix A{t) can be estimated roughly as 

Ajiit) r-u ^ t ^ oo 


and from (^51) it follows that at times of the order (^) the components 
Aji{t) with j,l ^ 1 are small compared with Pe^ and can be neglected in the 
diffusive exponent (|2^) . Thus, only An effectively survives. 

Let us then introduce a separation time to such that 1 (fD to -C InPe. 
In the limit Pe ^ oo this separation time will disappear from the hnal answer. 
The contribution to Q from the interval 0 < t < to is parametrically small, 
so that we can substitute 0 ^ to as a lower integration limit in 




Observe now that for the trajectories contributing to (Q”) the ^ij helds, 


being multiplied in the action by a factor exp(2pj — 2pj) 3> 1, will be 
exponentially depressed for f 3> 1. The components of T will therefore be 
frozen to constant values for those times. This is in particular true for the 
components rjj of the hrst row. 

Taking into account the freezing of T and writing explicitly A in the 
parametrization (p!8|, p!^ it is seen that the integral in ([TI|) saturates for f ~ to, 
so that the upper integration limit in ([TlD can be substituted with to- As a 
hnal result, in the limit Pe —> cx), (pH) asymptotically reduces to 




dt 


to 


d'^fc 

(27r)^ 




(27) 


X exp 




- 1 / 


(to)>V ^(r)>V ^’'^(r)T^'(to)no 


Observe that all the dependence of Q on the stochastic variables labelled by 
small times t < to is now concentrated in the square-bracketed exponent of 
(pTl) . The action (^) can also be written as a sum + 5^ of a small- 
and a large-times part. This makes clear the statement that integration over 
small-times variables, which describe a mainly dissipative process, can be 
separated from integration over large-times (t > to) variables, which describe 
multiplicative stretching. Time separation appears therefore as a natural 
consequence of the dynamic formalism. 

The small-times and large-times variables are coupled only by the condi¬ 
tion of continuity at t = to: 


Pj (^o) = pf{to), vf (to) = vf{to), X>^{to) = X<^{to). (28) 

In order to fully implement time separation we perform the reparametrization 

>V(t) ^>V(t)t-i(to). 

which implies the substitution T(to) —> 1 in (pT]) and modihes the boundary 
conditions (^^ID and (|^) as follows: 

p<(0) = 0, p<{to)=X<^{to) = 0, (29) 

for the small-times helds, and 


pj (to) = pf{to), vj (to) = X>^{to) = 0 (30) 
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for the large times. 

Thus, all the dissipative effects are encoded in the p.d.f. of the only 
variable 



(31) 


which appears as a parameter in the large-times averaging. In order to com¬ 
pute the total generating function 



(32) 


we shall integrate first on the large-times variables (essentially, only pi), find 
a p.d.f. which depends on /i as a parameter, and then complete the 
computation by integrating over the small-times p.d.f. This proce¬ 

dure can be suggestively thought as a renormalization group transformation 
consisting of only two steps. 

Let us now outline the first step. First of all, the variables can 

be integrated out completely, since they do not appear in (^). Of the 
variables, the only pf appears in (^7]). Integration over pf can be separated 
from the trivial integration over p^,..., p^_^ by the simple shift 



which also substitutes the initial condition pfito) = pfito) with the simpler 


Pj(to)=0, j = 


(34) 


Finally, we are left with an integration over the effective measure already 
obtained in Ref. 



(35) 


and with the factor 



( 36 ) 


which will be part of the small-times averaging measure. 
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4 Large-times averaging 


Let us now apply the previous considerations to the actual evaluation of 
dehned as the average of with respect to the measure (|35D with 
the constraint pi(to) = 0. The evaluation of the corresponding path integral 
can be reduced via the Feynman-Kac formula to the solution of the following 
auxiliary quantum mechanical problem in the space of the variable p = pi\ 


V: 


with 


lim e ^ 8 ^ *o) exp[—(T — 

T^oo 

lim ^(T - to; P = 0), (37) 

T^oo 


fe 

V>{p) 


(d-l)D 02 


dp^ 


sz(3 


D 


=2p 


+ ^>(P), 





/3 = , 


(38) 

(39) 

(40) 


In (^), the “wave function” T is the solution of the initial value problem 

= -dt^, (41) 

^(t = 0;p) = e^P. (42) 


For t —> cx) the behaviour of \h(t, p) is dominated by the lowest eigenvalue of 
the corresponding stationary problem, which can be easily determined since 
the potential V^{p) vanishes at p —>■ cx), leaving us with a free Hamiltonian: 


Ti^e^P ~ 


d {d — 1)2) d 


62'", p —i> CX) . 


So, for t —>• CX): 


d2(d-l)S) ^ d 

4/ (t; p) ~ 6 8 * e^P <F(p), 


(43) 
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with $(?/) the solution of the following boundary problem in the variable 
y = T)fiy/Pe^ e^: 


1 d rf+i 5 2sz „ 




0, (44) 


*^*( 1 /) ^ 1 , y ^ oo; y^{y)^0, y^O. 


(45) 


Note that <h(j/) can depend on the parameters s, z, X) and y only through 
their combination 


2s z 

{d — 1)2)/X ’ 


(46) 


and does not depend on Pe. 

Finally, from (|37|,|43|) it follows 


K. = $ 



(47) 


Here depends on Pe only through the argument of $. In the limit 
Fe —>■ cx) one simply hnds 

K. = < 1 >( 0 ) ( 48 ) 

For a generic given one can solve (^) iteratively in terms of the small 
potential sS(|/), s —0, hnding explicit expressions for all the “large-times” 
moments of Q\ 


K. = 1 + X;(2*:-1)!! 


k=l 


2sf3 


{d — 1)2)/i 


X 


(49) 


X 


/ k 

dx/i ■ ■ ■ dy2k n 

y2k < y2k+l 


y2i 

y2i-i 


d+l 


'Ay2i) 


The moments are here expressed directly in terms of the integral transform 
(^0[) of the pumping function Xk- 


5 Small times averaging 

In this section we shall compute 

Vf = {€-“■). (50) 
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The explicit expression of /i is 


pto 

/i = / V'^dt 

Jo 


with 

V = -h (51) 

In what follows we will omit the (.. .)^ index on the pj and pj variables. 
Along the lines previously exposed in the case of large-times averaging, the 
computation of Vf is again reduced to the solution of an auxiliary quantum 
mechanical problem: 


= exp 


- 1)D 


d-l 


x( H 6{X„ 


to X 


exp 


j=i 


2<m<n<d 


(52) 

d-l 

towo| !!'*(*))■ 

j=i 


Here 

n< = Ho+ sV, 

and Tdo is the quantum Hamiltonian corresponding to the classical action 

(H: 


Ho 


2 






dpidpj 


-dD ' 

l<J<2<d 


p2pi 2pj 


0<i<j<d ' ^ 


k<j; k<i 


dX, 


ki 


(53) 

(54) 


Due to the the triangular structure of the X^n and pj variables there arise 
no ordering problems. 

The seemingly intractable quantum mechanical problem in ^d{d+ 1) — 1 
dimensions can in fact be reduced to a solvable, one-dimensional problem 
thanks to the presence of |d(d -|- 1) — 2 symmetries. First of all, the initial 
wave function in (^) does not depend on the Xji variables. The evolution 
operator exp(—introduces a dependence on the pj variables, but not 
on the remaining Xij {i > 1), as is clear from the explicit expressions (^) 
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and (PH). We are thus led to consider the reduction of on the space 

of functions depending only on pj and pj\ 


d-l 






0<i<j<d 

i-1 


dpidpj 




d2 


dvLi 


(55) 


+ 1 /. (56) 


j=2 k=2 O- 

We can then forget altogether about the Xji variables and just compute 




'^5{pi) ■ ■■5{pd-i) exp ( - 6{pi) ■ ■ ■5(pd_i) 


(57) 


The reduced Hamiltonian depends on 2d —2 variables. Let us use now 

the fact that the action = 5^ + sp is invariant under the global left 
transformations 

T{t) ^ 0t(t), (58) 

where we take 0 = 1 + 50, with 50 having non zero values only in the hrst 
column and on the diagonal, and satisfying Tr(50) = 0 and (50)oo = 0. This 
way we generate 2d — 3 symmetries of S^, which are readily identihed with 
rotations in the space of the variables and shifts of the pj accompa¬ 

nied by corresponding rescalings of the pj. The corresponding variations of 
D and T can be read out from 

= (^5DD“^ - Dt"^5tD"^) (59) 

where (.. .)s denotes the symmetric part. The corresponding integrals of 
the motion are also easily found. They locally generate 2d — 3 “angular” 
directions, transverse to the “radial” direction labelled by V. We will not 
need their explicit expression: it will be enough to determine the “radial” 
part of 'H^P’P') by the transformation dx = dV/dx ■ dy + (radial part), with x 
substituted by pj, pj. In the to 3> 1 limit (|^ can be computed explicitly as 

/ oo poo 

dpi--- dpd-i ^pj = 0, Pj ), (60) 

-OO J —oo 
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p —>• CX) 


with <h the solution of the stationary problem 

= 0, 

d-l 

$(p,p) ~ 5(pi)---5(prf_i) exp j)pj 

i=i 

<h(p, p) ~ 0 , p —i> —oo . 

Let us consider a maximally symmetric solution <h having the form 

d-l 

$(p,p) =exp [^(d-i)pj 
j=i 


(61) 
(62) 

(63) 

(64) 

(65) 

(<i-l)S'- <“> 

where iLrf /2 is the Macdonald function of order d/2. The properties of the 


with T{y) satisfying the “radial” equation 


si/ 


2(d - + (d - l)(d + 2)VJ^' - —T = 0 . 


A solution of 


is: 


Tiy) = 


2sV 


d/2 


Macdonald functions |]^ allow to verify directly that (pif) satishes the bound¬ 
ary conditions (|52| , |53D . This also shows a posteriori that the boundary condi¬ 
tions (|6^j63|) are asymptotically invariant under the symmetries of the prob¬ 
lem. 

The integration (|60D can be explicitely performed, giving 


Vf = const ■ exp 


2s 


(d- 1)2) 

whose Laplace transform gives the p.d.f. of p: 

1 / I 


-P'-M = 


exp 


(67) 


( 68 ) 


y 2 ir(d-l)Sfi 3 /i V 2 (ci-l)B;i, 

One can check the soundness of the whole computation by verifying on the 
hrst moment of e that the constraint of energy conservation 


(^) = - 
' ' 2 


{2tiY 


■Xk 


is satished. 
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6 Resulting distribution function 


In order to restore the total p.d.f. V (e) we should take the function 
expressed in (^81) through the solution of then average over fi with the 
weight (^) and over z with the measure (|T^), and hnally perform an inverse 
Laplace transform. As already noticed, $ (?/ = 0) depends on the parameters 
s, z, D and fi only through their combination A (see (^6[)): 

$(0) = /(A) (69) 


The resulting V (e) can be expressed via / (A) as: 


V(e) 


27rr (d/2) y/e 


+ 00 



—oo 


da; 


+ 00 

y d^; / {ix) (-ix) exp 

0 




(70) 

The function / (A) has poles on the real negative semiaxis; the asymptotics 
of V {e) is dehned by the pole closest to the origin, which we denote as —A* 
(A* > 0): 


V (e) ~ exp 




oo 


(71) 


The point is that when A approaches the value —A* the Hamiltonian 
develops a ground state energy 


En < 


(E{d-l)'D 


and /(A) = given by (|37|) becomes inhnite. We therefore have 

“ 4 ^ 


with the maximum of the function y‘^'E{y) [cfr. (|38|,|39|) ). This hnally gives 


V (e) ~ exp —const (d^De) 


1/3' 


CXO . 


(72) 
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7 Discussion 


An important point to be noticed in the previous computation, expecially 
with regard to the asymptotic expression ([7^, is that we consider in this 
paper Pe as the largest number in the theory. For any large, but hnite Pe this 


means that the asymptotics (|7^ is valid ioi e / Eq ^ 1 (where £o = {^)) but 
smaller than InPe. There is here a crucial difference with the computation 
of the p.d.f. P{0) of the passive scalar 9 itself, which was discussed in the 
papers |^, ^ |^, [^. If 6* InPe the p.d.f. P{6) has a Gaussian form and 
can be obtained by taking into account only the variable pi corresponding to 
the highest Lyapunov exponent. However, if 0 > InPe, this approach does 
not give the correct exponent for d > 2, as it was observed in m- The 
correct asymptotic of P{9) for 6* S> InPe was computed in [|T^; the complete 
p.d.f was then restored in |T^, up to logarithmic corrections. 

Let us estimate the domain of validity of our expressions (^,^), which 
were obtained in the limit Pe —cx). In the computation, when pi ~ InPe 
we neglected the dependence of Q on p 2 ,...,pd-i- However, due to the 
peculiar dependence of the potential on pi, there exist trajectories for which 
Pi remains confined in a finite domain around InPe, while the variables pj, 
j = 2,... ,d — 1 freely evolve, eventually reaching values of the order In Pe. 
This happens for times of the order t* ~ In Pe. The corresponding value of 
e is eiim ~ Eot* ~ eolnPe. For e > eum the behaviour of V{£) is probably 
modified. 

The exponent of the p.d.f.’s tail (^) up to d-dependent numerical coef¬ 
ficient coincides with the one-dimensional results |^. This means that the 
configurations of the passive scalar field contributing into (|7^) are effectively 
one-dimensional, as it was noted in . 
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A Computation of the Jacobian 

In this appendix we derive the form of the Jacobian term (|^) . Let us 
introduce a basis ejj of the matrix algebra gl((J) dehned by 


make use of the notation 

tp j ^ b dfc Gfcfc) k 1, . . . d, fjj Qij Gji) 3 ^ ^ 1 

and parametrize 

, D = exp (^pfcdfcj , = exp {XijUj) , 

where summation on repeated indices is assumed and Pd = — (pi + - ■ ■ + Pd-i)- 
Under a variation 6pk, SXij) the matrix a = VVVV“^ varies as 

Jd = A + [A, a] = crlfiij + cr^'^^dfc + alfUj, with A = dW ■ 

The measure Va is invariant under the global symmetries W —> TZqW and 
VV —>■ so that the Jacobian 


R = T exp 


J 


det 


5a 


J(0,p,X) 


does not depend on the R and T variables. Setting then R = 1, T = 1 after 
having computed the variation 5a we hnd 

+ id~^5(j)ij){pi - pj) , 

(^if = -ef'^~P^5Xij-2{d~^5(j)ij){pi-pj). 

The Jacobian matrix has thus triangular form and can be computed with the 
help of the regularization 
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( 73 ) 





The result (^) is finally obtained making use of the simple identity 

d d—1 

~ Pi^ = + ^^Pj = 2 “ ^^Pi ■ 

l<i<j<d j=l j=l 

Note that the consistency of the regularization prescription (0) can be 
checked by computing ([>V(f)no]^) first directly and that using 
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